Improving polygenic risk prediction in admixed populations by explicitly modeling ancestral-differential effects via GAUDI.

Polygenic risk scores (PRS) have shown successes in clinics, but most PRS methods focus only on participants with distinct primary continental ancestry without accommodating recently-admixed individuals with mosaic continental ancestry backgrounds for different segments of their genomes. Here, we develop GAUDI, a novel penalized-regression-based method specifically designed for admixed individuals. GAUDI explicitly models ancestry-differential effects while borrowing information across segments with shared ancestry in admixed genomes. We demonstrate marked advantages of GAUDI over other methods through comprehensive simulation and real data analyses for traits with associated variants exhibiting ancestral-differential effects. Leveraging data from the Women's Health Initiative study, we show that GAUDI improves PRS prediction of white blood cell count and C-reactive protein in African Americans by > 64% compared to alternative methods, and even outperforms PRS-CSx with large European GWAS for some scenarios. We believe GAUDI will be a valuable tool to mitigate disparities in PRS performance in admixed individuals.

Analytical challenges in omics research on asthma and allergy: A National Institute of Allergy and Infectious Diseases workshop.

Studies of asthma and allergy are generating increasing volumes of omics data for analysis and interpretation. The National Institute of Allergy and Infectious Diseases (NIAID) assembled a workshop comprising investigators studying asthma and allergic diseases using omics approaches, omics investigators from outside the field, and NIAID medical and scientific officers to discuss the following areas in asthma and allergy research: genomics, epigenomics, transcriptomics, microbiomics, metabolomics, proteomics, lipidomics, integrative omics, systems biology, and causal inference. Current states of the art, present challenges, novel and emerging strategies, and priorities for progress were presented and discussed for each area. This workshop report summarizes the major points and conclusions from this NIAID workshop. As a group, the investigators underscored the imperatives for rigorous analytic frameworks, integration of different omics data types, cross-disciplinary interaction, strategies for overcoming current limitations, and the overarching goal to improve scientific understanding and care of asthma and allergic diseases.

Session Introduction: Overcoming health disparities in precision medicine.

The following sections are included:OverviewDealing with the lack of diversity in current research datasetsDevelopment of fair machine learning algorithmsRace, genetic ancestry, and population structureConclusionAcknowledgments.

A statistical framework for powerful multi-trait rare variant analysis in large-scale whole-genome sequencing studies.

Large-scale whole-genome sequencing (WGS) studies have improved our understanding of the contributions of coding and noncoding rare variants to complex human traits. Leveraging association effect sizes across multiple traits in WGS rare variant association analysis can improve statistical power over single-trait analysis, and also detect pleiotropic genes and regions. Existing multi-trait methods have limited ability to perform rare variant analysis of large-scale WGS data. We propose MultiSTAAR, a statistical framework and computationally-scalable analytical pipeline for functionally-informed multi-trait rare variant analysis in large-scale WGS studies. MultiSTAAR accounts for relatedness, population structure and correlation among phenotypes by jointly analyzing multiple traits, and further empowers rare variant association analysis by incorporating multiple functional annotations. We applied MultiSTAAR to jointly analyze three lipid traits (low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides) in 61,861 multi-ethnic samples from the Trans-Omics for Precision Medicine (TOPMed) Program. We discovered new associations with lipid traits missed by single-trait analysis, including rare variants within an enhancer of NIPSNAP3A and an intergenic region on chromosome 1.

Mosaic chromosomal alterations in blood across ancestries using whole-genome sequencing.

Megabase-scale mosaic chromosomal alterations (mCAs) in blood are prognostic markers for a host of human diseases. Here, to gain a better understanding of mCA rates in genetically diverse populations, we analyzed whole-genome sequencing data from 67,390 individuals from the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine program. We observed higher sensitivity with whole-genome sequencing data, compared with array-based data, in uncovering mCAs at low mutant cell fractions and found that individuals of European ancestry have the highest rates of autosomal mCAs and the lowest rates of chromosome X mCAs, compared with individuals of African or Hispanic ancestry. Although further studies in diverse populations will be needed to replicate our findings, we report three loci associated with loss of chromosome X, associations between autosomal mCAs and rare variants in DCPS, ADM17, PPP1R16B and TET2 and ancestry-specific variants in ATM and MPL with mCAs in cis.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed whole-genome sequencing study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs.

Late-life plasma proteins associated with prevalent and incident frailty: A proteomic analysis.

Proteomic approaches have unique advantages in the identification of biological pathways that influence physical frailty, a multifactorial geriatric syndrome predictive of adverse health outcomes in older adults. To date, proteomic studies of frailty are scarce, and few evaluated prefrailty as a separate state or examined predictors of incident frailty. Using plasma proteins measured by 4955 SOMAmers in the Atherosclerosis Risk in Community study, we identified 134 and 179 proteins cross-sectionally associated with prefrailty and frailty, respectively, after Bonferroni correction (p < 1 × 10-5 ) among 3838 older adults aged ≥65 years, adjusting for demographic and physiologic factors and chronic diseases. Among them, 23 (17%) and 82 (46%) were replicated in the Cardiovascular Health Study using the same models (FDR p < 0.05). Notably, higher odds of prefrailty and frailty were observed with higher levels of growth differentiation factor 15 (GDF15; pprefrailty = 1 × 10-15 , pfrailty = 2 × 10-19 ), transgelin (TAGLN; pprefrailty = 2 × 10-12 , pfrailty = 6 × 10-22 ), and insulin-like growth factor-binding protein 2 (IGFBP2; pprefrailty = 5 × 10-15 , pfrailty = 1 × 10-15 ) and with a lower level of growth hormone receptor (GHR, pprefrailty = 3 × 10-16 , pfrailty = 2 × 10-18 ). Longitudinally, we identified 4 proteins associated with incident frailty (p < 1 × 10-5 ). Higher levels of triggering receptor expressed on myeloid cells 1 (TREM1), TAGLN, and heart and adipocyte fatty-acid binding proteins predicted incident frailty. Differentially regulated proteins were enriched in pathways and upstream regulators related to lipid metabolism, angiogenesis, inflammation, and cell senescence. Our findings provide a set of plasma proteins and biological mechanisms that were dysregulated in both the prodromal and the clinical stage of frailty, offering new insights into frailty etiology and targets for intervention.

Evaluating genomic signatures of aging in brain tissue as it relates to Alzheimer's disease.

Telomere length (TL) attrition, epigenetic age acceleration, and mitochondrial DNA copy number (mtDNAcn) decline are established hallmarks of aging. Each has been individually associated with Alzheimer's dementia, cognitive function, and pathologic Alzheimer's disease (AD). Epigenetic age and mtDNAcn have been studied in brain tissue directly but prior work on TL in brain is limited to small sample sizes and most studies have examined leukocyte TL. Importantly, TL, epigenetic age clocks, and mtDNAcn have not been studied jointly in brain tissue from an AD cohort. We examined dorsolateral prefrontal cortex (DLPFC) tissue from N = 367 participants of the Religious Orders Study (ROS) or the Rush Memory and Aging Project (MAP). TL and mtDNAcn were estimated from whole genome sequencing (WGS) data and cortical clock age was computed on 347 CpG sites. We examined dementia, MCI, and level of and change in cognition, pathologic AD, and three quantitative AD traits, as well as measures of other neurodegenerative diseases and cerebrovascular diseases (CVD). We previously showed that mtDNAcn from DLPFC brain tissue was associated with clinical and pathologic features of AD. Here, we show that those associations are independent of TL. We found TL to be associated with ß-amyloid levels (beta = - 0.15, p = 0.023), hippocampal sclerosis (OR = 0.56, p = 0.0015) and cerebral atherosclerosis (OR = 1.44, p = 0.0007). We found strong associations between mtDNAcn and clinical measures of AD. The strongest associations with pathologic measures of AD were with cortical clock and there were associations of mtDNAcn with global AD pathology and tau tangles. Of the other pathologic traits, mtDNAcn was associated with hippocampal sclerosis, macroscopic infarctions and CAA and cortical clock was associated with Lewy bodies. Multi-modal age acceleration, accelerated aging on both mtDNAcn and cortical clock, had greater effect size than a single measure alone. These findings highlight for the first time that age acceleration determined on multiple genomic measures, mtDNAcn and cortical clock may have a larger effect on AD/AD related disorders (ADRD) pathogenesis than single measures.

Whole Genome Sequencing Based Analysis of Inflammation Biomarkers in the Trans-Omics for Precision Medicine (TOPMed) Consortium.

Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38,465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program. We identified 22 distinct single-variant associations across 6 traits - E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin - that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed Whole Genome Sequencing Study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions. Large-scale whole genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess the associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with blood lipid levels (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare variant aggregate association tests using the STAAR (variant-Set Test for Association using Annotation infoRmation) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare coding variants in nearby protein coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500 kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variations and rare protein coding variations at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNA, implicating new therapeutic opportunities.

Age at Menopause, Leukocyte Telomere Length, and Coronary Artery Disease in Postmenopausal Women.

BACKGROUND: Premature menopause is a risk factor for accelerated cardiovascular aging, but underlying mechanisms remain incompletely understood. This study investigated the role of leukocyte telomere length (LTL), a marker of cellular aging and genomic instability, in the association of premature menopause with cardiovascular disease. METHODS: Participants from the UK Biobank and Women's Health Initiative with complete reproductive history and LTL measurements were included. Primary analyses tested the association between age at menopause and LTL using multivariable-adjusted linear regression. Secondary analyses stratified women by history of gynecologic surgery. Mendelian randomization was used to infer causal relationships between LTL and age at natural menopause. Multivariable-adjusted Cox regression and mediation analyses tested the joint associations of premature menopause and LTL with incident coronary artery disease. RESULTS: This study included 130 254 postmenopausal women (UK Biobank: n=122 224; Women's Health Initiative: n=8030), of whom 4809 (3.7%) had experienced menopause before age 40. Earlier menopause was associated with shorter LTL (meta-analyzed ß=-0.02 SD/5 years of earlier menopause [95% CI, -0.02 to -0.01]; P=7.2×10-12). This association was stronger and significant in both cohorts for women with natural/spontaneous menopause (meta-analyzed ß=-0.04 SD/5 years of earlier menopause [95% CI, -0.04 to -0.03]; P<2.2×10-16) and was independent of hormone therapy use. Mendelian randomization supported a causal association of shorter genetically predicted LTL with earlier age at natural menopause. LTL and age at menopause were independently associated with incident coronary artery disease, and mediation analyses indicated small but significant mediation effects of LTL in the association of menopausal age with coronary artery disease. CONCLUSIONS: Earlier age at menopause is associated with shorter LTL, especially among women with natural menopause. Accelerated telomere shortening may contribute to the heightened cardiovascular risk associated with premature menopause.

Incorporating genetics in identifying peanut allergy risk and tailoring allergen immunotherapy: A perspective on the genetic findings from the LEAP trial.

Examining the genetics of peanut allergy (PA) in the context of clinical trial interventions and outcomes provides an opportunity to not only understand gene-environment interactions for PA risk but to also understand the benefit of allergen immunotherapy. A consistent theme in the genetics of food allergy is that in keeping with the dual allergen exposure hypothesis, barrier- and immune-related genes are most commonly implicated in food allergy and tolerance. With a focus on PA, we review how genetic risk factors across 3 genes (FLG, MALT1, and HLA-DQA1) have helped delineate distinct allergic characteristics and outcomes in the context of environmental interventions in the Learning Early about Peanut Allergy (LEAP) study and other clinical trials. We specifically consider and present a framework for genetic risk prediction for the development of PA and discuss how genetics, age, and oral consumption intertwine to predict PA outcome. Although there is some promise in this proposed framework, a better understanding of the mechanistic pathways by which PA develops and persists is needed to develop targeted therapeutics for established disease. Only by understanding the mechanisms by which PA develops, persists, and resolves can we identify adjuvants to oral immunotherapy to make older children and adults immunologically similar to their younger, more malleable counterparts and thus more likely to achieve long-term tolerance.

A functional genomics pipeline to identify high-value asthma and allergy CpGs in the human methylome.

BACKGROUND: DNA methylation of cytosines at cytosine-phosphate-guanine (CpG) dinucleotides (CpGs) is a widespread epigenetic mark, but genome-wide variation has been relatively unexplored due to the limited representation of variable CpGs on commercial high-throughput arrays. OBJECTIVES: To explore this hidden portion of the epigenome, this study combined whole-genome bisulfite sequencing with in silico evidence of gene regulatory regions to design a custom array of high-value CpGs. This study focused on airway epithelial cells from children with and without allergic asthma because these cells mediate the effects of inhaled microbes, pollution, and allergens on asthma and allergic disease risk. METHODS: This study identified differentially methylated regions from whole-genome bisulfite sequencing in nasal epithelial cell DNA from a total of 39 children with and without allergic asthma of both European and African ancestries. This study selected CpGs from differentially methylated regions, previous allergy or asthma epigenome-wide association studies (EWAS), or genome-wide association study loci, and overlapped them with functional annotations for inclusion on a custom Asthma&Allergy array. This study used both the custom and EPIC arrays to perform EWAS of allergic sensitization (AS) in nasal epithelial cell DNA from children in the URECA (Urban Environment and Childhood Asthma) birth cohort and using the custom array in the INSPIRE [Infant Susceptibility to Pulmonary Infections and Asthma Following RSV Exposure] birth cohort. Each CpG on the arrays was assigned to its nearest gene and its promotor capture Hi-C interacting gene and performed expression quantitative trait methylation (eQTM) studies for both sets of genes. RESULTS: Custom array CpGs were enriched for intermediate methylation levels compared to EPIC CpGs. Intermediate methylation CpGs were further enriched among those associated with AS and for eQTMs on both arrays. CONCLUSIONS: This study revealed signature features of high-value CpGs and evidence for epigenetic regulation of genes at AS EWAS loci that are robust to race/ethnicity, ascertainment, age, and geography.

Structural variation across 138,134 samples in the TOPMed consortium.

Ever larger Structural Variant (SV) catalogs highlighting the diversity within and between populations help researchers better understand the links between SVs and disease. The identification of SVs from DNA sequence data is non-trivial and requires a balance between comprehensiveness and precision. Here we present a catalog of 355,667 SVs (59.34% novel) across autosomes and the X chromosome (50bp+) from 138,134 individuals in the diverse TOPMed consortium. We describe our methodologies for SV inference resulting in high variant quality and >90% allele concordance compared to long-read de-novo assemblies of well-characterized control samples. We demonstrate utility through significant associations between SVs and important various cardio-metabolic and hematologic traits. We have identified 690 SV hotspots and deserts and those that potentially impact the regulation of medically relevant genes. This catalog characterizes SVs across multiple populations and will serve as a valuable tool to understand the impact of SV on disease development and progression.

Structural variation across 138,134 samples in the TOPMed consortium.

Ever larger Structural Variant (SV) catalogs highlighting the diversity within and between populations help researchers better understand the links between SVs and disease. The identification of SVs from DNA sequence data is non-trivial and requires a balance between comprehensiveness and precision. Here we present a catalog of 355,667 SVs (59.34% novel) across autosomes and the X chromosome (50bp+) from 138,134 individuals in the diverse TOPMed consortium. We describe our methodologies for SV inference resulting in high variant quality and >90% allele concordance compared to long-read de-novo assemblies of well-characterized control samples. We demonstrate utility through significant associations between SVs and important various cardio-metabolic and hemotologic traits. We have identified 690 SV hotspots and deserts and those that potentially impact the regulation of medically relevant genes. This catalog characterizes SVs across multiple populations and will serve as a valuable tool to understand the impact of SV on disease development and progression.

The IL-4Rα Q576R polymorphism is associated with increased severity of atopic dermatitis and exaggerates allergic skin inflammation in mice.

BACKGROUND: Atopic dermatitis (AD) is characterized by TH2-dominated skin inflammation and systemic response to cutaneously encountered antigens. The TH2 cytokines IL-4 and IL-13 play a critical role in the pathogenesis of AD. The Q576->R576 polymorphism in the IL-4 receptor alpha (IL-4Rα) chain common to IL-4 and IL-13 receptors alters IL-4 signaling and is associated with asthma severity. OBJECTIVE: We sought to investigate whether the IL-4Rα R576 polymorphism is associated with AD severity and exaggerates allergic skin inflammation in mice. METHODS: Nighttime itching interfering with sleep, Rajka-Langeland, and Eczema Area and Severity Index scores were used to assess AD severity. Allergic skin inflammation following epicutaneous sensitization of mice 1 or 2 IL-4Rα R576 alleles (QR and RR) and IL-4Rα Q576 (QQ) controls was assessed by flow cytometric analysis of cells and quantitative RT-PCR analysis of cytokines in skin. RESULTS: The frequency of nighttime itching in 190 asthmatic inner-city children with AD, as well as Rajka-Langeland and Eczema Area and Severity Index scores in 1116 White patients with AD enrolled in the Atopic Dermatitis Research Network, was higher in subjects with the IL-4Rα R576 polymorphism compared with those without, with statistical significance for the Rajka-Langeland score. Following epicutaneous sensitization of mice with ovalbumin or house dust mite, skin infiltration by CD4+ cells and eosinophils, cutaneous expression of Il4 and Il13, transepidermal water loss, antigen-specific IgE antibody levels, and IL-13 secretion by antigen-stimulated splenocytes were significantly higher in RR and QR mice compared with QQ controls. Bone marrow radiation chimeras demonstrated that both hematopoietic cells and stromal cells contribute to the mutants' exaggerated allergic skin inflammation. CONCLUSIONS: The IL-4Rα R576 polymorphism predisposes to more severe AD and increases allergic skin inflammation in mice.

MUC5B, telomere length and longitudinal quantitative interstitial lung changes: the MESA Lung Study.

BACKGROUND: The MUC5B promoter variant (rs35705950) and telomere length are linked to pulmonary fibrosis and CT-based qualitative assessments of interstitial abnormalities, but their associations with longitudinal quantitative changes of the lung interstitium among community-dwelling adults are unknown. METHODS: We used data from participants in the Multi-Ethnic Study of Atherosclerosis with high-attenuation areas (HAAs, Examinations 1-6 (2000-2018)) and MUC5B genotype (n=4552) and telomere length (n=4488) assessments. HAA was defined as the per cent of imaged lung with attenuation of -600 to -250 Hounsfield units. We used linear mixed-effects models to examine associations of MUC5B risk allele (T) and telomere length with longitudinal changes in HAAs. Joint models were used to examine associations of longitudinal changes in HAAs with death and interstitial lung disease (ILD). RESULTS: The MUC5B risk allele (T) was associated with an absolute change in HAAs of 2.60% (95% CI 0.36% to 4.86%) per 10 years overall. This association was stronger among those with a telomere length below an age-adjusted percentile of 5% (p value for interaction=0.008). A 1% increase in HAAs per year was associated with 7% increase in mortality risk (rate ratio (RR)=1.07, 95% CI 1.02 to 1.12) for overall death and 34% increase in ILD (RR=1.34, 95% CI 1.20 to 1.50). Longer baseline telomere length was cross-sectionally associated with less HAAs from baseline scans, but not with longitudinal changes in HAAs. CONCLUSIONS: Longitudinal increases in HAAs were associated with the MUC5B risk allele and a higher risk of death and ILD.

Recent progress in the genetic and epigenetic underpinnings of atopy.

In the past 2 years, there continue to be advances in our understanding of the genetic and epigenetic underpinnings of atopy pertaining to disease risk and disease severity. The joint role of genetics and the environment has been emphasized in multiple studies. Combining genetics with family history, biomarkers, and comorbidities is further refining our ability to predict the development of individual atopic diseases as well as the advancement of the atopic march. Polygenic risk scores will be an important next step for the field moving toward clinical translation of the genetic findings thus far. A systems biology approach, as illustrated by studies of the microbiome and epigenome, will be necessary to fully understand disease development and to develop increasingly targeted therapeutics.

Session Introduction: Overcoming health disparities in precision medicine.

The following sections are included: Overview, Equitable risk prediction, Pharmacoequity, Race, genetic ancestry, and population structure, Conclusion, Acknowledgments, References.